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1.0  INTRODUCTION 


The  interaction  between  the  number  of  snapshots  (temporal  sa.np'es)  and  the 
number  of  sensors  in  adaptive  beamforming  is  recognized  in  literature.  It  is  well  known 
from  the  seminal  paper  of  Reed,  Mallett,  and  Brennan  (1974)  that,  compared  to  the 
ideal  case  of  an  infinite  number  of  snapshots,  the  loss  in  signal-to-interference-plus- 
noise  ratio  of  the  Minimum  Variance  Distortionless  Response  (MVDR)  is  less  than  3 
dB  if  the  number  of  snapshots  of  data  is  at  least  twice  the  number  of  array  elements. 
More  fundamentally,  when  the  number  of  snapshots  is  fewer  than  the  number  of 
sensors,  the  sample  cross-spectral  density  matrix  R  (which  is  also  the  maximum 
likelihood  estimator  of  the  true,  but  unknown,  cross-spectral  density  matrix  R)  is 
singular,  rendering  the  MVDR  method  useless.  This  led  Capon  (1969)  to  modify 
the  cross-spectral  density  matrix  R  by  adding  a  small  amount  of  incoherent  noise  to 
the  diagonal  of  R;  the  procedure  now  is  known  as  diagonal  loeuling  or  white-noise 
stabilization  (Carlson,  1988;  Tran,  1990). 

For  a  very  large  array  of  size  n  and  a  relatively  small  number  of  snapshots  T, 
the  modified  MVDR  method  encounters  another  problem,  namely  the  computational 
load  in  inverting  the  (n  x  n)  modified  cross-spectral  density  matrix.  The  recently 
introduced  JASON  algorithm  (Rothaus,  1991)  requires  the  inversion  of  a  (T  x  T) 
matrix,  a  desirable  advantage  especially  when  T  <^n. 

The  purpose  of  this  report  is  to  investigate  properties  and  performance  of  conven¬ 
tional  (c  oartlett),  MVDR  (with  diagonal  loading),  and  JASON  beamformers  when 
the  numbi !  of  snapshots  is  small  compared  to  the  number  of  array  elements. 

The  remainder  of  the  report  is  orga’Mzed  v.to  five  ''•eotions.  In  secti  --n  2.0,  a  nar¬ 
rowband,  multiple-source  model  is  presented.  The  formulation  is  an  extension  of  a 
narrowband,  single-sour- e  model  presented  in  Rothaus  (1991).  This  section  provides 
notations  and  terminology  for  section  3.0  where  a  quick  review  of  the  three  beamform¬ 
ing  methods,  conventional,  M’\'DR,  and  JASON  beamformers,  is  given  with  attention 
to  the  situation  where  the  number  of  snapshots  is  small  relative  to  the  number  of  ar¬ 
ray  elements.  In  section  4.0,  tii-’  general  model  in  section  2.0  is  refined  for  simulation. 
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In  particular,  we  consider  an  environment  consisting  of  twenty  narrowband,  plane- 
wave  sources,  corrupted  in  additive  colored  noise,  impinging  on  a  uniform  line  array 
of  100  elements  with  interelement  spacing  equal  to  one  half  of  the  array  wavelength. 
The  noise  is  zero  mean  and  may  be  temporally  correlated  from  sample  to  sample. 
In  section  5.0,  comparative  simulation  results  of  the  three  beamforming  methods  are 
presented  as  a  function  of  the  number  of  snapshots  that  is  less  than  or  equal  to  the 
number  of  sensors.  Qualitative  results  based  on  output  array  response,  and  quantita¬ 
tive  results  b2ised  on  normalized  output  deflection  and  peak-to-background  ratio,  are 
presented  in  subsections  5.1  and  5  2,  respectively.  Finally,  in  section  6.0,  we  conclude 
the  report  by  sumarizing  our  findings  with  emphasis  on  a  well-populated  array  and 
a  small  number  of  snapshots. 

An  interactive  code,  written  in  MATLAB,  for  plotting  array  response  is  included 
in  Appendix  A. 

2.0  GENERAL  MODEL 

In  this  section  a  general  model  is  formulated  for  a  narrowband,  multisource  en¬ 
vironment  after  Rothaus  (1991)  where  a  single  source  is  explicitly  considered.  Most 
notations  and  terminology  in  Rothaus  (1991)  are  followed  for  quick  reference. 

Suppose  we  have  an  acoustic  field  of  n  sensors  distributed  in  space,  a  point  source 
in  the  direction  of  the  steering  vector  a  and  k  other  sources,  called  interferences, 
with  direction  specified  by  steering  vectors  fij,  j  =  I,'-  ■  ,k. 

At  time  t  and  sensor  t,  by  changing  the  local  time  at  each  sensor,  the  observed 
data  can  be  written  as 

Xi{t)  =  NM  +  OiPit)  +  (1) 

}=i 

where  Ni{i)  denotes  the  noise  field,  P{t)  and  Qj{t),  j  =  1,- ••,!•,  are  acoustic  pres¬ 
sures,  or  complex  magnitudes  (envelopes),  of  sources. 

By  stacking  expression(l)  vertically,  serially  in  i,  for  t  =  1,2,. . .  ,n,  we  arrive  at 
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the  following  compact  matrix  notation 


X{t)  =  N(t)  +  Pit)a  +  (2) 

3=1 

The  vector  X(t)  is  often  referred  to  as  the  snapshot  at  time  t. 

Consider  T  consecutive  snapshots  X(<),  X{t  +  1),  X{t  +  T  —  1).  Assume 
that  the  steering  vectors  a,  j  =  1,. . . ,  A:,  remain  constant  during  the  processing 
interval  represented  by  the  duration  of  the  T  snapshots.  Let 


S{i)  ^  [x(0:X(t  +  l):  •••  iX(t  +  7’-l)], 
Afit)  t  [N(oiN(<  +  i) ...  :N(<  +  r-i)], 


Pit) 
Pit  + 1) 


^(0 


i  = 


+  r  - 1) 


TTiit  +  T-l) 


Then,  we  obtain 

Sit)  =  Afit)  +  aVit)  +  ZPiQ-it).  (3) 

i=i 

Here,  the  “overbar”  denotes  complex  conjugate,  and  *,  conjugate  transpose. 

It  follows  immediately  that  R,  the  maximum-likelihood  estimate  of  the  true  (but 
unknown)  cross-spectral  density  matrix  Rof  the  received  pressure  field  can  be  written 
as 

—  A  1  1 

R  =  f  E  x(’-)x-(t)  =  yS(i)s-(<).  (4) 

Notice  that  the  average  signal  power  is  ^  V*it)Vit)  and,  similarly,  the  average 
power  of  the  j-th  interference  is  ^  Qjit)Qjit).  Also,  the  average  noise  power  per 
sensor  is 

where  tr(Y)  is  the  trace  of  matrix  Y. 

The  space  spanned  by  columns  of  S,  which  are  the  T  snapshots,  is  called  the 
observation  space.  These  snapshots  can  be  assumed  to  be  linearly  independent. 
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3.0  BEAMFORMING 


A  quick  review  of  three  beamforming  methods,  conventiorai,  MVDR,  and  the 
recently  introduced  JASON,  is  presented  in  this  section.  Our  emphasis  is  on  the 
situation  where  the  number  of  snapshots  is  relatively  small  compared  to  the  number 
of  sensors. 

The  MVDR  beamformer  is  an  adaptive,  high- resolution  beamformer  that  mini¬ 
mizes  the  output  power  while  maintaining  unity  response  in  the  direction  of  look. 
Mathematically,  a  weight  vector  w  is  sought  after  for  the  constrained  optimization 
problem 

min  w*Rw  subject  to  w*7  =1,  (5) 


where  7  denotes  the  look-direction  steering  vector. 
The  optimal  weight  vector  is  readily  obtained  as 

R-*7 
7*R-*7’ 


w  = 


which  yields  the  optimal  output  power 

w*Rw  = 


1 


7.R-I7 

The  conventional  beamformer,  described  by  the  delay-and-sum  process,  is  equiv¬ 
alent  to  setting  the  weight  vector  equal  to  the  look-direction  steering  vector  7.  Its 
output  power  is  7*R7. 

Operationally,  due  to  the  unknown  nature  of  the  cross-spectral  density  matrix  R, 
we  deal  exclusively  with  its  maximum-likelihood  estimator  R.  Thus,  by  replacing  R 
by  R  and  using  equation  (4),  an  operational  form  for  problem  (5)  is 


niin  w*55*w  subject  to  w*7  =  1.  (6) 

Here,  the  constant  term  ^  is  ignored  since  it  has  no  effect  on  the  minimization. 

By  taking  expectation,  we  obtain  the  “true"  output  powers  corresponding  to  T 
snapshots: 

POWERmvob  =  e(^)  =  (7) 
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and 


POWERo.  =  E(7-R7)  =  (7-55-7)  ■  (8) 

Since  the  n  x  n  matrix 

_  I  t+T-l 

R  =  i  E  X(’-)X‘W 

■*  T=l 

is  a  scalar  multiple  of  the  sum  of  T  matrices  each  having  rank  one,  the  rank  of  R  is  at 
most  T.  Therefore,  when  the  number  of  snapshots  T  is  smaller  than  n,  the  number  of 
sensors,  R  is  necessarily  singular,  thus  rendering  equation  (7)  useless  (Capon,  1969). 

In  practice,  however,  when  T  <  n,  white-noise  stabilization  (diagonal  loading)  is 
applied  to  obtain  the  invertibility  for  R  (Capon,  1%9).  The  method  calls  for  adding 
to  R  (prior  to  inverting)  a  diagonal  matrix  that  can  be  written  in  the  form 


where  tr  is  the  trace  operator,  I  the  identity  matrix,  and  the  stabilization  factor  F 
is  typically  chosen  between  0.001  to  0.1  (Tran,  1990).  Notice  that  a  stabilization 
factor  0.001  corresponds  to  adding  to  the  system  a  white  sensor  noise  30  dB  below 
the  average  sensor  power  level. 

The  JASON  algorithm  (Rothaus,  1991)  offers  an  alternative  to  the  white-noise 
stabilization  method  in  solving  the  constrained  optimization  problem  (5),  or  (6), 
when  the  number  of  snapshots  is  less  than  the  number  of  sensors.  The  algorithm  is 
based  on  the  observation  that  if  a  point  source  is  present  in  the  direction  specified  by 
7,  and  if  the  noise-plus-interference  vector 

M(r)  4  N(r)  +  (9) 

i=i 

in  the  snapshot  X(t),  for  r  =  +  T  —  1,  has  all  mean  values  zero,  then  the 

desired  weight  vector  w  lies  in  the  range  space  of  E  (5*5) .  To  show  this,  we  let 
7  =  0,  where  a  is  the  steering  vector  of  a  point  source  with  acoustic  pressures 
F(0, P{t  +  T-l).  Let 

M(()  4  ^-(0  + 

j=I 
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Then 


w*E(55*)w  =  w*E(>(M*)w  +  ^{V^V)  |w•a|^ 

Since  E(A<A^*)  is  positive  semidefinite,  and  since  the  scalar  ^{V*V)  >  0  (sig¬ 
nal  present),  the  assumption  that  w  belongs  to  the  null  space  of  E  {S*S) ,  that  is 
E  {S*S)  w  =  0,  implies  immediately  w*a  =  0.,  Therefore,  to  satisfy  problem  (5)  with 
7  =  0;,  the  weight  vector  w  must  lie  in  the  range  space  of  E  {SS”) . 

Operationally,  the  range  space  of  E  (55*)  can  be  taken  to  be  the  space  spanned 
by  the  columns  of  5,  that  is  the  observation  space.  Therefore,  the  JASON  algorithm 
solves  the  constrained  optimization  problem  (6)  with  an  additional  requirement 

w  belongs  to  the  observation  space.  (10) 

Consequently,  we  can  represent  the  weight  vector  w  as  a  linear  combination  of  T 
snapshots.  That  is, 

w  =  5c, 

where  c  is  the  vector  of  size  T  whose  entries  are  linear  combination  coefficients. 
Thus,  instead  of  solving  problem  (6)  for  the  weight  vector  w,  one  solves  the  following 
transformed  problem  for  the  coefficient  vector  c 


min  c*5*55*5c  subject  to  c*5*7  =  1. 


(11) 


The  optimal  solution  is 


(5*5)-' 5*7 


c  = 


7*5  (5*5)-^  5*7’ 

which  yields  the  estimated  output  power  corresponding  to  T  snapshots 

w*Rw  =  4 ^ — 5 - .  (12) 

T  7*5  (5*5)"' 5*7  ^  ’ 

The  “true”  output  power  of  the  JASON  algorithm  corresponding  to  T  snapshots  is 

given  by 


POWERjason  =  is 


1 


f  -  Vr5(5*5)-^^j  • 

Notice  that  the  T  x  T  matrix  5*5  is  invertible  since  the  T  snapshots  are  linearly 
independent. 
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The  advantage  of  the  JASON  algorithm  lies  in  the  invertibility  of  S*S,  in  contrast 
to  the  singularity  of  the  matrix  when  the  number  of  snapshots  T  is  less  than 
the  number  of  sensors  n.  Furthermore,  when  T  n,  as  in  the  case  of  a  populated 
array,  inverting  the  T  xT  matrix  S*S  is  less  a  computational  burden  than  inverting 
any  n  x  n  matrix  obtained  from  5<S*  via  diagonal  modification. 

4.0  SIMULATION  MODEL 

In  this  section,  the  general  model  in  section  2.0  is  refined  for  simulation.  In  par¬ 
ticular,  we  consider  an  environment  consisting  of  20  narrowband,  plane-wave  sources, 
corrupted  by  additive  noise,  impinging  on  a  uniform  line  array  of  100  elements  with 
interelement  spacing  equal  to  one-half  of  the  wavelength.  Operational  details  of  the 
simulation  are  also  presented  in  this  section. 

For  the  snapshot, 

X(0  =  N(0  -b  Pit)a  -f  (14) 

j=i 

we  wish  to  model  noise  vector  N(t),  direction  vectors  a,  j  =  and  the 

corresponding  complex  magnitudes  P{t),  Qj{t),  j  =  1, 2, . . . ,  fc. 

Such  modeling  will  be  applied  to  the  T  consecutive  snapshots  X(i),  X(t  -f  1),  . . . , 
X(t  -b  r  -  1),  as  in 

S(t)  =  m  +  aV{t)  +  (15) 

J  =  1 

where  we  recall 

S{t)  ^  [x(t)  i X(f -b  1) i  •••  iX(t-br-i)], 

^  [n(o  :  N(f  -b  1)  •  •  •  i  N(t  -b  r  - 1)1 ,  (16) 

and 

m  ]  ^(0 

+  1)  ^  A  Q}{i  +  1) 

P{i  +  T-1)  J  + 
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We  shall  consider  the  case  where  sources  located  in  the  farfield  can  be  modeled 
as  plane  waves.  Also,  for  simplicity  we  shall  consider  a  line  array  of  n  elements, 
uniformly  distributed  with  interelement  spacing  equal  to  one-half  of  the  array  wave¬ 
length. 

For  a  line  array,  the  angle-of-arrival  of  a  plane  wave  is  mezisured  relative  to  the 
normal  to  the  array.  Consequently,  the  range  of  the  angle-of-arrival  is  —90°  to  90°. 

We  recall  the  following  useful  result  (Haykin,  1985).  For  the  coordinate  reference 
at  the  center  of  the  array,  the  direction  vector  of  a  plane  wave  impinging  on  the 
uniform  linear  array  of  n  elements  at  an  angle-of-arrival  <^,  measured  with  respect  to 
the  normal  to  the  array,  is  given  by 


. , , , , 
1,  ^ . 

where  superscript  T  stands  for  transpose,  i  =  and 

27rd  . 

w  =  — r-  sm(9) 

A 


n  even, 
n  odd. 


with  d  as  the  interelement  spacing  and  A  the  array  wavelength.  Note  that  the  “center” 
of  the  array  is  its  actual  middle  sensor  if  n  is  odd,  otherwise  it  is  a  fictitious  center. 
Therefore,  with  the  coordinate  reference  at  the  center  of  the  array,  the  direction 
vector  of  a  plane  wave  impinging  on  the  uniform  linear  array  can  be  obtained  readily 
by  setting  d  =  A/2. 

In  the  following  we  will  realize  W(f),  V{t)y  Qj{t),  j  =  l,...,fc,  using  normally 
distributed  random  numbers  available  in  the  MATLAB  package. 

The  function  rand(n,T)  in  MATLAB  generates  an  (n  x  r)-matrix  with  random 
entries.  By  specifying  the  function  rand  as  normal  beforehand,  the  entries  of  the  re¬ 
sulting  matrix  will  be  normally  distributed  with  mean  0.0  and  variance  1.0.  Similarly, 


rand(n,T)  +  i  rand(n,T)  (17) 

generates  a  complex  random  matrix  of  size  nxT  whose  real  part  and  imaginary  part 
of  each  of  the  entries  are  normally  distributed  with  means  0.0  and  variances  1.0.  This 
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can  be  taken  as  a  prototype  for  M{t)  and,  similarly, 

rand(T,  1)  +  i  rand(T,  l)  (18) 

for  'P(t),  and  Qj{t),  j  =  1, . . . ,  Observe  that  the  complex  vector  in  expression  (18) 
and  columns  of  the  complex  matrix  in  (17)  are  realizations  of  zero  mean,  complex 
Gaussian  random  vectors. 

For  simulation  we  wish  to  generate  M{t),  V{t),  Qj{t),  j  =  1,...,A:,  as  random 
matrices  with  zero  means  and  with  average  noise  power  per  sensor  <t^,  average  signal 
power  ffg,  and  average  power  of  the  j-th.  interference  cr|,  where  we  recall, 

<^1  4  itr(iAf(i)Af(())  =  itrWi)V(i)), 

"I  =  f'P'imt),  <>]  =  j  =  \x-x 

Furthermore,  for  simplicity  we  choose  <7^  =  1 . 

For  any  matrix  A  =  (a,j) ,  recall  that 

tr{AA-)  =  =  IIAfp, 

hi 

where  ||A||p  denotes  the  Frobenius  norm  of  A.  In  particular,  we  can  write 

Consider  the  normalization 

(rand(n,T)  +  i  rand(n,T)) 
||rand(n,T)Firand(n,T)llp  ^  ^ 

Then  the  average  noise  power  is  unity. 

Strictly  speaking,  the  noise  is  no  longer  Gaussian  since  the  denominator  of  equa¬ 
tion  (19)  is  randomly  distributed.  The  distribution  of  ||rcuid(n,T)  -f  i  rand(n,T)||p 
depends  largely  on  the  statistical  independence  (or  dependence)  of  the  2nT  normally 
distributed  random  variables  (each  with  mean  0.0  and  variance  1.0)  that  make  up 
the  real  and  imaginary  parts  of  entries  of  the  matrix  rand(n,T)  +  i  rand(n,T).  In  the 
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ideal  situation  where  they  are  statistically  independent,  then  the  denominator  of  (1 9) 
has  the  chi  distribution  (Melsa  &  Sage,  1973,  pp.  314-315). 

Without  precise  knowledge  of  the  distribution  of  ||rand(n,T) -f  i  raiid(n,T)|jp ; 
however,  for  practical  purposes,  we  may  assume  its  expectation  exists  and  is  greater 
than  zero.  The  assurr  ption  on  the  positivity  of  the  expected  value  of  the  denom¬ 
inator  of  equation  (19)  is  justified  since  ||rand(n, T) -f  i  raiid(n, T)||p  represents  a 
magnitude- type”  quantity  as  signified  by  the  definition  of  Frobenius  norm.  Con¬ 
sequently, 

p,  ^ (rand(n,T)  +  i  ra]id(n,T))\  E  (y/nT  (rand(n,T)  +  i  rEmd(n,T))) 

I  ||rand(n,T)-Hirand(n,T)||p  )  E(|lrand(n,T)  -|- i  rand(n,T)||p)  ^ 
that  is,  M{t)  is  a  complex  random  matrix  with  mean  0„xr- 

In  view  of  equation  (16),  the  noise  is  zero  mean.  Furthermore,  it  is  colored  since 
it  may  be  temporally  correlated  from  sample  to  sample. 

Similarly,  we  consider 

V(t)  =  (rand(T,l)  -f  irand(T,l)) 

||rand(T,  1)  -|-  irand(T,  1)||  ’  ^  ' 

where  ll'H  stands  for  the  usual  Euclidean  norm,  a  special  case  of  the  Frobenius  norm. 
Then  ^  P‘‘(t)P(t)  =  cr].  Moreover,  it  follows  from  the  above  argument  that  V(t)  is 
a  complex  random  vector  with  mean  Otxi- 

Finally,  Qj(i),  j  =  can  be  generated  by 

^  /air  (rand(T,l)  -f-  i  rand(T,  1)) 

^  ||rand(T,l)  -h  i  rand(T,l)||  ’ 

As  a  consequence  of  the  modeling  of  equations  (19),  (20)  and  (21),  for  any  source 
chosen  as  the  desired  signal,  the  noise-plus-interference  vector  M(r)  =  N(t)  + 
Sj=i  Qj{'^)(^j  of  the  data  snapshot  X(t ),  for  t  =  <,...,/  -t-  T  —  1 ,  has  all  mean  values 
of  zero. 


For  illustration,  an  interactive  MATLAB  implementation  for  plotting  array  re¬ 
sponse  is  presented  in  Appendix  A.  This  code  can  be  run  on  a  PC.  Actual  simulations, 
however,  were  done  on  a  CONVEX  C220,  a  two-processor  minisupercomputer  at  the 
Naval  Ocean  Systems  Center. 
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5.0  SIMULATION  RESULTS 


In  this  section,  we  describe  the  multisource  environment  used  in  simulation  to¬ 
gether  with  qualitative  and  quantitative  results  concerning  properties  and  perfor¬ 
mance  of  conventional,  MVDR  (with  stabilization  factors  0.001,  0.01,  and  0.1),  and 
JASON  beamformers  when  the  number  of  snapshots  is  small.  Qualitative  results, 
based  on  output  array  response,  and  quantitative  results,  based  on  normalized  out¬ 
put  deflection  and  peak-to-background  ratio,  are  presented  in  subsections  5.1  and  5.2, 
respectively. 

We  consider  20  sources,  corrupted  by  additive  noise,  impinging  on  a  line  array  of 
100  equally  spaced  sensors  with  interelement  spacing  equal  to  one  half  of  the  array 
wavelength.  The  noise  components  N(Q,...,N(f  -I-  T  -  1)  in  the  T  consecutive 
snapshots  are  modeled  by  equation  (19)  as  being  random  with  zero  means  and  unit 
average  power  per  sensor.  Direction  vectors  of  sources  are  assumed  to  remain  constant 
during  the  processing  interval  indicated  by  the  duration  of  the  T  snapshots.  Acoustic 
pressures,  or  complex  amplitudes,  of  sources  are  modeled  as  in  equation  (20)  where 
average  source  powers  are  specified  relative  to  average  noise  power  per  sensor  (0  dB). 
Table  1  below  lists  source  angles  and  the  corresponding  powers. 


Table  1.  Angle  and  power  of  20  sources. 


source 

angle  (deg) 

source  power  (dB) 

1 

-65 

-5 

2 

-45 

-15 

3 

-30 

-10 

4 

-10 

5 

5-  15 

10,  12,  14,  16, 

18,  20,  22,  24, 

26,  28,  30 

0 

16-20 

40,  42,  44, 

46  48 
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Notice  that  there  are  two  clusters  of  sources  of  0-dB  average  power;  the  first  group 
is  comprised  of  11  sources,  the  second  of  5  sources,  with  2°  spacing  in  each  group. 
This  is  the  smallest  angular  separation  between  any  2  of  the  20  sources. 

The  Rayleigh  limit  of  angular  resolution  of  a  uniform  linear  array  of  n  elements  is 
approximately  2/(n  — 1)  radians  (Krolik  &:  Swingler,  1989).  In  particular,  for  n  =  100 
sensors,  the  Rayleigh  limit  is  approximately  1.1575°,  which  is  less  than  the  2°  spacing 
just  mentioned. 

5.1  Qualitative  Results 

The  array  response  of  the  conventional  beamformer,  of  MVDR  with  stabilization 
factors  0.001,  0.01,  0.1,  and  of  the  JASON  algorithm  are  plotted  for  the  look-direction 
steering  vector  sweeping  from  —90°  to  90°  with  1°  increments.  The  number  of 
snapshots  considered  are  10,  20,  30,  40  and  50,  that  is,  from  to  |  of  the  number 
of  elements. 

Corresponding  to  each  number  of  snapshots,  we  plotted  the  array  response  of 
five  independent  trials  and,  in  a  separate  graph,  the  average  response  taken  over 
200  independent  trials  for  each  method.  These  plots  of  the  array  response  provide  a 
qualitative  assessment  of  the  properties  and  performance  of  each  beamformer. 

For  a  quick  reference,  figures  of  array  response  of  algorithms  are  listed  in  tables 
2  through  6.  For  example,  table  2  shows  that  figure  l(al)  depicts  conventional  array 
response  of  5  independent  trials  whereeis  figure  l(a2)  exhibits  the  average  response  of 
the  conventional  beamformer  taken  over  200  independent  trials,  and  that  these  two 
figures  were  obtained  in  the  case  where  the  number  of  snapshots  T  =  10. 
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Table  2.  List  of  figures  of  conventional  array  response. 


number  of  snapshots 

10 

20 

30 

40 

50 

Figure 

l(al) 

five  ind 

l(bl) 

ependei 

l(cl) 

it  trials 

l(dl) 

l(el) 

ge  of  2C 

i(b2) 

0  indep 

l(c2) 

endent 

l(d2) 

trials 

l(e2) 

Table  3.  List  of  figures  of  MVDR  array  response, 
(stabilization  factor  0.001). 


number  of  snapshots 

10 

20 

30 

40 

50 

five  independent  trials 

Figure 

2(al) 

2{cl) 

2(dl) 

2(el) 

average  of  200  independent  trials 

2(a2) 

2(b2) 

2(c2) 

2(d2) 

2(e2) 

Table  4.  List  of  figures  of  MVDR  array  response, 
(stabilization  factor  0.01). 


number  of  snapshots 

10 

20 

30 

40 

50 

five  independent  trials 

Figure 

3(al) 

3{bl) 

3(cl) 

3(dl) 

3(el) 

average  of  200  independent  trials 

3(a2) 

3(b2) 

3(c2) 

3(d2) 

3{e2) 

Table  5.  List  of  figures  of  MVDR  array  response, 
(stabilization  factor  0.1). 


Table  6.  List  of  figures  of  JASON  array  response. 


number  of  snapshots 

10 

20 

30 

40 

50 

five  independent  trials 

Figure 

5(al) 

5(bl) 

5(cl) 

5(dl) 

5(el) 

average  of  200  independent  trials 

5(a2) 

6{b2) 

5(c2) 

5(d2) 

5(e2) 
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-90  -60  -30  0  30  60  90 

angle  (deg) 


(al)  Five  independent  trials,  each  with  10  snapshots. 


(a2)  Average  of  200  independent  trials,  each  with  10  snapshots 


Figure  1.  Conventional  array  response. 
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angle  (deg) 


(b2)  Average  of  200  independent  trials,  each  with  20  snapshots. 
Figure  1  (cont.).  Conventional  array  response. 
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(d2)  Average  of  200  independent  trials,  each  with  40  snapshots. 


Figure  1  (cont.).  Conventional  array  response. 
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angle  (deg) 

(el)  Five  independent  trials,  each  with  50  snapshots. 


elevation  (deg) 

(e2)  Average  of  200  independent  trials,  each  with  50  snapshots. 
Figure  1  (cont.).  Conventional  array  response. 
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(bl)  Five  independent  trials,  each  with  20  snapshots. 


Figure  2  (cent.).  MVDR  array  response  with  stabilization  factor  0.001. 
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(cl)  Five  independent  trials,  each  with  30  snapshots. 


(c2)  Average  of  200  independent  trials,  each  with  30  snapshots. 


Figure  2  (cont.).  MVDR  array  response  with  stabilization  factor  0.001 


(dl)  Five  independent  trials,  each  with  40  snapshots. 


angle  (deg) 

(d2)  Average  of  200  independent  trials,  each  with  40  snapshots. 


Figure  2  (cont.).  MVDR  array  response  with  stabilization  factor  0.001 


Figure  2  (coiit.).  MVDR  array  response  with  stabilization  factor  0.001. 
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(al)  Five  independent  trials,  each  with  10  snapshots. 


Figure  3.  MVDR  array  response  with  stabilization  factor  0.01. 
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/j 


Figure  3  (cont.).  MVDR  array  response  with  stabilization  factor  0.01. 
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angle  (deg) 


(cl)  Five  independent  trials,  each  with  30  snapshots. 


angle  (deg) 


(c2)  Average  of  200  independent  trials,  each  with  30  snapshots. 
Figure  3  (cont.).  MVDR  array  response  with  stabilization  factor  0.01. 
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angle  (deg) 


(dl)  Five  independent  trials,  each  with  40  snapshots. 


angle  (deg) 


(d2)  Average  of  200  independent  trials,  each  with  40  snapshots. 
Figure  3  (cont.).  MVDR  array  response  with  stabilization  factor  0.01. 
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(e2)  Average  of  200  independent  trials,  each  with  50  snapshots. 


Figure  3  (cont.).  MVDR  array  response  with  stabilization  factor  0.01. 
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(al)  Fiye  independent  trials,  each  with  10  snapshots. 


(a2)  Average  of  200  independent  trials,  each  with  10  snapshots. 
Figure  4.  MVDR  array  response  with  stabilization  factor  0.1. 
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• 

' 

i 

j 

i 

. 

1 

: 

i 

angle  (deg) 

(bl)  Five  independent  trials,  each  with  20  snapshots. 


ULA  olilOO  sensors 


angle  (deg) 

(b2)  Average  of  200  independent  trials,  each  with  20  snapshots. 


Figure  4  (cont,).  MVDR  array  response  with  stabilization  factor  O.j 


(cl)  Five  independent  trials,  each  with  30  snapshots. 


(c2)  Average  of  200  independent  trials,  e2w:h  with  30  snapshots. 


Figure  4  (cont.).  MVDR  array  response  with  stabilization  factor  0.1. 
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Figure  4  (coat.).  MVDR  array  response  with  stabilization  factor  0.1. 
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-90  -60  -30  0  30  60  90 

angle  (deg) 


(al)  Five  independent  trials,  each  with  10  snapshots. 


(a2)  Average  of  200  independent  trials,  each  with  10  snapshots. 


Figure  5.  JASON  array  response. 
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•90  >60  -30  0  30  60  90 

angle  (deg) 

(bl)  Five  independent  trials,  each  with  20  snapshots. 


-90  -60  -30  0  30  60  90 

angle  (deg) 

(b2)  Average  of  200  independent  trials,  each  with  20  snapshots. 


Figure  5  (cont.).  JASON  array  response. 


(cl)  Five  independent  trials,  each  with  30  snapshots. 


(c2)  Average  of  200  independent  trials,  each  with  30  snapshots. 


Figure  5  (cont.).  JASON  array  response. 
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(dl)  Five  independent  trials,  each  with  40  snapshots. 


Figure  5  (cont.).  JASON  array  response. 
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5 


(e2)  Average  of  200  independer*  trials,  each  with  50  snapshots. 


Figure  5  (cont.).  JASON  array  response. 
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Peaks  in  these  figures  of  array  response  should  be  interpreted  in  relative  terms 
with  respect  to  the  corresponding  background  level.  For  example,  for  the  5-dB  source 
at  —10°  with  40  snapshots,  the  200-trial  average  output  power  relative  to  background 
level  is  about  24  dB  for  the  conventional  beamformer  (figure  l(d2)),  about  19.5  dB  for 
the  MVDR  with  stabilization  factor  0.001  (figure  2(d2)),  about  19  dB  for  the  MVDR 
with  stabilization  factor  0.01  (figure  3(d2)),  about  16.25  dB  for  the  MVDR  with 
stabilization  factor  0.1  (figure  4(d2)),  and  about  13.5  dB  for  the  JASON  beamformer 
(figure  4(d2)). 

The  overall  observation  is  that,  for  each  method,  the  higher  the  number  of  snap¬ 
shots,  the  lower  the  variability  of  estimate.  Moreover,  the  higher  number  of  snapshots 
also  reduces  the  variability  of  background  level  relative  to  the  primary  peaks.  For 
the  MVDR  beamformer,  the  higher  the  stabilization  factor,  the  lower  the  primary 
peaks  relative  to  background  level,  and  the  higher  the  stabilization  factor,  the  lower 
the  variability  of  background  level  relative  to  the  primary  peaks.  Finally,  among  the 
three  methods,  the  JASON  beamformer  tends  to  have  the  highest  level  of  variability 
both  in  estimates  and  in  background  level  (relative  to  primary  peaks). 

For  each  value  of  T  =  10,20,30,40,50,  the  conventional  beamformer  tends  to 
outperform  other  methods  as  it  resolves  each  source  with  the  highest  peak  relative 
to  background  level.  This  is  not  surprising  since  sources  are  well  spaced  with  the 
smallest  angular  separation  greater  than  the  Rayleigh  limit.  Indeed,  for  a  very  large 
uniform  array,  the  Rayleigh  angular  limit  is  significantly  reduced,  hence  for  all  prac¬ 
tical  purposes,  resolution  is  no  longer  working  again  t  the  conventional  beamformer 
eis  in  the  case  of  a  small  uniform  array. 

When  the  number  of  snapshots  is  small,  we  observe  an  inverted  behavior  of  the 
JASON  array  response,  namely,  yielding  a  dip  instead  of  a  peak  at  the  source’s 
location.  This  inverted  behavior  is  dependent  on  powers  of  sources.  As  the  number 
of  snapshots  increases,  sources  with  high  SNRs  begin  to  peak  in  the  usual  manner 
earlier  than  those  with  low  SNRs. 
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To  understand  the  inverted  behavior  of  the  JASON  array  response  when  the  num¬ 
ber  of  snapshots  is  small,  let  us  consider  the  extreme  ctise  of  one  snapshot.  For 
r  =  1,  then  the  observation  matrix  S  reduces  to  a  vector,  therefore  S*S  =  ||«5||^ 
is  a  (positive)  scalar.  Consequently,  the  output  power  of  the  JASON  beamformer 
corresponding  to  R,  a  realization  of  R,  is  given  by 

w*Rw  =  (22) 

7*R7 

Notice  that  the  denominator  of  equation  (22)  is  the  estimated  power  of  the  conven¬ 
tional  beamformer  corresponding  to  R.  When  the  look- direction  steering  vector  7 
coincides  with  the  steering  vector  of  any  source,  7*R7  peaks  and,  hence,  its  recipro¬ 
cal  yields  a  dip  at  the  arrival  angle.  Clearly  the  numerator  {S*S)^  of  equation  (22) 
is  independent  of  the  look-direction  steering  vector. 

The  foregoing  discussion,  indeed,  can  be  recast  by  a  short,  but 
nonconstructive,  argument:  for  T  =  1, 

POWER--  = 

\  7*R7  /  E  (7*R7)  POWERcon 
The  inverted  behavior  of  the  JASON  array  response  when  T  is  small,  therefore, 
can  be  explained  as  an  extrapolation  of  the  case  of  one  snapshot. 

5.2  Quantitative  Results 

To  quantify  output  array  response  of  each  beamformer,  we  consider  the  normalized 
output  deflection  (Burdic,  1978;  Krolik  &  Swingler,  1990)  that  is  defined  by 

DEFLECTION  (POWER(^)) 

where  POWER(^)  stands  for  estimated  output  power  when  steered  at  the  angle 
(f)  and  Std(D)  denotes  standard  deviation  of  □.  Clearly,  a  high-normalized  output 
deflection  is  a  desirable  property  for  detection. 


E(POWER(?i)  I  signal  present^  —  E  (POWER(<^)  |  signal  absent) 
Std  ^POWER(^)  I  signal  absent) 


,(23) 
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To  assess  the  relative  detection  performance  of  each  beamformer  when  the  number 
of  snapshots  is  small,  with  interest  in  the  detectability  of  a  weak  source  in  the  presence 
of  strong  interferences,  we  compute  normalized  output  deflection  of  each  algorithm 
when  the  look-direction  elevation  angle  is  steered  at  —45°,  the  direction  of  the  weakest 
(—15  dB)  among  the  twenty  sources,  for  T  running  from  10  to  100  with  increment 
equal  to  5  snapshots.  Each  deflection  estimate  is  obtained  from  200  independent 
trials,  using  sample  means  of  output  power  when  steered  at  —45°  either  signal  is 
present  or  not,  and  standard  deviation  of  output  power  when  steered  at  —45°  in  the 
absence  of  signal. 

For  real  situations  of  multiple  sources,  the  normalized  output  deflection  defined  as 
in  equation  (23)  is  difficult  to  estimate  since  one  cannot  recreate  the  same  environment 
except  for  one  source  that  now  is  absent.  Furthermore,  equation  (23)  is  evaluated 
at  one  particular  location  and  ignores  the  background  that  might  contain  spurious 
peaks. 

An  alternative  to  the  normalized  output  deflection  is  the  peak-to-background  ra¬ 
tio,  defined  as 

peak  output  power  —  average  background  power  . 

standard  deviation  of  the  baickground  power  ' 

The  mean  background  output  level  is  usually  calculated  excluding  a  small  interval 
around  the  peak  (Feuillade,  DelBalzo,  &  Rowe,  1989).  Obviously,  for  a  power  detec¬ 
tor,  higher  peak-to-background  ratio  indicates  better  detection. 

Other  things  being  equal,  we  note  that,  from  equation  (24),  higher  variability  of 
background  level  results  in  lower  peak-to-background  ratio  and,  from  equation  (23), 
higher  variability  of  estimates  results  in  lower  normalized  output  deflection. 

For  the  environment  of  multiple  sources  depicted  in  table  1,  we  compute  the  peak- 
to-background  ratio  for  the  look-direction  angle  —45°  as  follows.  The  average  and 
standard  deviation  of  background  level  is  calculated  over  the  18  grid  points  from 
-55°  to  -35°  with  1°  spacing,  excluding  output  power  at  -46°,  -45°  and  -44°. 
Peak-to-background  ratio  for  each  beamformer  when  steered  at  —45°  is  obtained  as 
the  ensemble  mean  of  estimates  from  200  independent  trials.  For  comparison,  we 
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consider  the  number  of  snapshots  from  10  to  100  with  increments  of  5  snapshots, 
that  is,  from  ^  of  the  number  of  sensors  to  the  number  of  sensors. 

Because  the  main  lobe  of  the  conventional  beamformer  is  relatively  wide  compared 
to  that  of  the  MVDR  and  the  JASON,  if  the  interval  excluded  from  computing  the 
mean  background  output  level  is  too  small,  there  might  be  a  downward  bias  in  the 
resulting  peak-to-background  ratio  for  the  conventional  beamformer.  The  effect  of 
the  length  of  the  excluded  interval,  however,  depends  on  the  number  of  elements.  For 
the  uniform  linear  array  n  elements,  one  measure  of  separation  is  the  first-null  (or 
standard)  beamwidth,  which  is  defined  as  the  angular  separation  between  the  peak 
of  the  main  lobe  and  the  first  null  of  the  radiation  pattern  function  lF(a;)  = 
where  u  =  ^sin(^)  with,  we  recall,  d  as  the  interelement  spacing,  A  the  array 
wavelength,  and  <f>  the  steering  angle.  Let  BW  denote  a  first-null  beamwidth,  then 
BW  =  ^  radians  or  ^  degrees  (Haykin,  1985).  In  our  ceise,  the  excluded  interval  is 
the  open  interval 

(-45°  -  A^,-45°  -f  A(^)  , 

where  A^  =  2°.  For  the  array  of  100  equally-spaced  elements  with  interelement 
spacing  equal  to  one-half  of  the  array  wavelength,  A<^  =  2°  amounts  to 
Aw  =  50  sin(7r/90)  BW  =  1.V45  BW,  and  hence  the  immediate  grid  points  outside 
the  excluded  interval,  —47°  and  —43°,  are  well  separated  from  the  main  lobe. 

The  normalized  output  deflection  and  peak-to-background  ratio  of  each  beam- 
former  is  plotted,  on  linear  scale  and  on  dB  scale,  versus  the  number  of  snapshots. 
Due  to  the  inverted  behavior  of  the  JASON  array  response  when  the  number  of 
snapshots  is  small  relative  to  the  number  of  elements,  we  consider  absolute  values 
of  normalized  output  deflection  and  peak-to-background  ratio  of  the  JASON  beam- 
former.  This  amounts  to  taking  absolute  values  of  the  numerators  of  equations  (23) 
and  (24)  in  computing  these  expressions.  This  provision  is  satisfactory  for  detection 
cissessment. 
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Overall  results  of  normalized  output  deflection  and  peak-to-background  ratio  are 
essentially  the  same.  The  only  noticeable  difference  between  these  two  metrics  occurs 
for  the  MVDR,  as  a  group,  when  the  number  of  snapshots  is  under  25.  Figure  6(b) 
shows  that  the  normalized  deflection  of  the  MVDR  for  each  of  the  three  different 
stabilization  factors  is  approximately  8  dB,  4  dB,  and  3  dB  below  the  normalized 
deflection  of  the  conventional  beamformer  for  T  =  10,  15,  and  20,  respectively.  For 
these  values  of  T,  figure  7(b)  indicates  that  the  peak-to-background  ratio  of  the  con¬ 
ventional  beamformer  is  less  than  2  dB  above  that  of  the  MVDR  for  each  stabilization 
factor.  This  discrepancy  might  be  attributed  to  the  fact  noted  earlier  that  the  normal¬ 
ized  output  deflection  is  computed  when  steered  at  a  particular  angle  (—45°)  without 
taking  into  consideration  the  existence  of  spurious  peaks  in  the  background  of  the 
source.  When  T  is  as  low  as  10,  15  or  20,  the  relatively  less  variable  background 
of  the  MVDR  beamformer  (figures  2(al)  to  2(b2),  figures  3(al)  to  3(b2),  and  figures 
4(al)  to  4(b2)),  accounts  for  small  standard  deviation  of  background  level,  and  thus 
higher  peak-to-background  ratio  compared  to  the  corresponding  normalized  output 
deflection. 

Results  in  figures  6(a)  and  7(a)  show  that,  for  the  MVDR  beamformer,  there  is  an 
indication  that  higher  stabilization  factor  yields  higher  deflection  and  higher  peak-to- 
background  ratio,  and  the  effect  increases  as  the  number  of  snapshots  approaches  the 
number  of  elements.  This  is  consistent  with  earlier  qualitative  observations  that  both 
higher  number  of  snapshots  and  higher  stabilization  factor  promote  smaller  standard 
deviation  of  estimate  and  smaller  standard  deviation  of  background  level  relative  to 
primary  peaks.  When  the  number  of  snapshots  T  <  50,  on  the  dB  scale,  figure 
6(b)  indicates  that,  for  the  MVDR  beamformer  with  stabilization  factors  0.001,  0.01, 
and  0.1,  the  differences  in  normalized  output  deflection  are  within  1  dB.  The  same 
observation  is  noted  for  the  MVDR’s  peak-to-background  ratios  in  figure  7(b)  when 
T  <  50. 
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number  of  snapshots 

Figure  6(a).  Normalized  output  deflection  (on  linear  scale, 
with  absolute  value  for  JASON)  versus  number  of  snapshots 
for  the  source  with  angle-of- arrival  —45°. 
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Figure  6(b).  Normalized  output  deflection  (on  linear  scale, 
with  absolute  value  for  JASON)  versus  number  of  snapshots 
for  the  source  with  angle-of-arrival  -45°. 
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Figure  7(a).  Peak-to-background  ratio  (on  linear  scale, 
with  absolute  value  for  JASON)  versus  number  of  snapshots 
for  the  .source  with  angle-of-arrival  -45°. 
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peak-tO“background  ratio  (dB) 


Figure  V^b).  Peak-to-background  ratio  (on  dB  scale, 
with  absolute  value  for  JASON)  versus  number  of  snapshots 
for  the  source  with  angle-of-arrival  —45°. 
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For  the  JASON  beamformer,  when  the  number  of  snapshots  is  small  compared  to 
the  number  of  sensors,  all  in  all,  the  low  normalized  output  deflection  is  the  result 
of  high  variability  of  the  estimate;  similarly  its  low  peak-to-background  ratio  is  a 
result  of  high  variability  of  background  level.  The  dip  at  T  =  45  in  figure  6(b) 
and  figure  7(b)  is  a  superficial  artifact  resulting  from  the  of  taking  absolute  values 
of  normalized  deflection  and  peak-to-background  ratio.  Indeed,  for  T  <  40,  both 
normalized  deflection  and  peak-to-background  ratio  are  negative  and  relatively  larger, 
in  absolute  values,  than  the  corresponding  one  at  T  =  45.  Therefore,  taking  absolute 
values  of  normalized  deflection  and  peak-to-background  ratio  yields  a  mirror  image 
of  each  when  T  <  40  and  a  dip  at  T  =  45,  as  in  figures  6(a)  and  7(a);  this  effect  is 
accentuated  on  the  dB  scale,  cis  in  figures  6(b)  and  7(b). 

6.  CONCLUSION 

The  properties  and  performance  of  conventional.  Minimum  Variance  Distortionless 
Response  (MVDR),  and  JASON  beamformers  were  investigated  via  simulation  for  a 
multiple-source  environment  corrupted  by  colored  noise  with  special  emphasis  on  the 
situation  where  the  number  of  snapshots  is  small  compared  to  the  number  of  array 
elements.  All  in  all,  for  a  uniform  linear  array  of  100  elements  with  interelement 
spacing  equal  to  one-half  of  the  array  wavelength,  when  sources  are  well  separated, 
the  conventional  beamformer  tends  to  outperform  the  other  two  methods  as  it  resolves 
each  source  with  the  highest  peak  relative  to  background  level.  With  attention  to 
large  arrays,  the  choice  of  well  separated  sources  in  the  study  is  justified.  Indeed,  for 
a  very  large  uniform  array,  the  Rayleigh  angular  limit  is  significantly  reduced.  Hence 
for  all  practical  purposes,  resolution  is  no  longer  working  against  the  conventional 
beamformer  as  in  the  case  of  a  small  uniform  array.  The  MVDR  beamformer,  with 
white-noise  stabilization,  displays  better  resolution  exemplified  by  sharp  peaks  at  the 
expense  of  inverting  of  the  modified  cross-spectral  density  matrix,  a  computational 
burden  for  a  well-populated  array.  For  the  MVDR  beamformer,  we  further  notice 
that  the  higher  the  stabilization  factor,  the  lower  the  variability  of  background  level 
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relative  to  the  primary  peaks  as  well  as  the  lower  the  primary  peaks  relative  to 
background  level.  The  recently  introduced  JASON  beamformer  exhibits  an  inverted 
behavior,  namely  yielding  a  dip  instead  of  a  peak  at  source  location,  when  the  number 
of  snapshots  is  small  compared  to  the  number  of  sensors.  This  inverted  behavior  of 
the  JASON  array  response  can  be  explained  as  an  extrapolation  of  the  behavior  for 
the  special  case  of  one  snapshot.  Moreover,  the  behavior  is  more  persistent  for  sources 
with  low  SNRs.  As  the  number  of  snapshots  increases,  sources  with  high  SNRs  begin 
to  peak  in  the  usual  manner  earlier  than  those  with  low  SNRs. 

Quantitative  results,  bzised  on  normalized  output  deflection  and  peak-to  background 
ratio,  indicated  better  detection  performance  of  conventional  and  MVDR  beamform 
ers.  When  the  number  of  snapshots  is  small  relative  to  the  number  of  array  elements, 
the  low  normalized  output  deflection  of  the  JASON  algorithm  is  caused  by  the  high 
variability  of  the  estimate,  and  similarly  its  low  peak-to-background  ratio  is  a  result 
of  the  high  variability  of  background  level. 

Findings  presented  in  this  report  were  based  on  the  simulation  using  an  array 
model  with  relatively  low-average  sidelobes.  Some  improvements  of  the  performance 
of  the  adaptive  algorithms,  especially  the  modified  MVDR,  relative  to  conventional 
beamforming  might  be  expected  for  an  array  model  with  high  sidelobes. 
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APPENDIX  A 

MATLAB  Code  for  Plotting  Array  Response 


A-1 


For  illustration,  an  interactive  MATLAB  code  for  plotting  array  response  of  a 
uniform  linear  array  is  presented  in  this  appendix. 

The  following  code  is  runnable  on  8086/8088-based  computers  (PC/XT)  with  PC- 
MATLAB,  or  on  80286-based  computers  (AT)  with  AT-MATLAB.  However,  these 
versions  limit  vectors  to  a  length  of  8188  elements  and  square  matrices  to  order  90. 
It  is  preferable  to  run  the  code  on  80386-based  computers  with  386- MATLAB  where 
variables  can  have  unlimited  size,  subject  to  only  the  ammount  of  memory  installed 
and  the  amount  of  free  hard-disk  space. 

Upon  execution,  the  code  will  ask  for  inputs  including  the  following: 

the  number  of  sensors  n, 

the  number  of  snapshots  T, 

the  number  of  (plane  wave)  sources  k, 

angle-of-arrival  (AOA)  and  signal-to-noise  ratio  (SNR)  of  each  source,  and 
the  number  of  independent  trials. 

The  user  is  able  to  choose  among  JASON,  MVDR,  and  conventional  beamformers 
for  plotting  array  response.  For  MVDR,  the  user  is  prompted  to  specify  a  value  for 
stabilization  factor  between  0.001  and  0.1.  Clearly,  for  T  >  n,  a  *’alue  of  zero  can  be 
supplied,  provided  the  resulting  cross-spectral  density  matrix  is  i  '.11-conditioned.. 

The  code  then  computes  and  plots  the  array  response  of  each  independent  trial  and 
the  corresponding  average,  together  on  the  same  axes  of  the  first  plot.  For  JASON 
algorithm,  due  to  its  high  variability  of  estimate,  peaks  of  subsequent  trials  may  be 
out  of  range  of  the  first  plot,  especially  when  n  is  small  and  T  is  smaller  than  n.  We 
did  not  attempt  to  adjust  the  axes  to  accommodate  all  trials. 
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Columns  of  the  matrix  below  are  steering  vectors  of  k  sources 


o 

oo 


a 

* 

•  ^ 

a 

.a 

‘53 


a 

* 


+ 

“O 

a 


(N 


CO 

a 

,2 

(H 

•73 

A! 

O 

O 


u 

0) 

I 

3 

3 


a. 

O 

o 


0 

3 

* 

1 

HH 

3 

00 

hJ 

3 

04 

r; 

rH 

D 

3 

i-H  ^ 

II 

II 

II  11 

CO 

CO 

».< 

4) 

X  ti 

& 

3  -3 

’2  O'  « 

a  •— <  V.O 

s 

?  w 

9  vQ  vO 

3 

3 

0 

a  .S 
«  3 

6  g 

3  3 

CO  PJ 

s 

■w  p 

s  « 

'C  k 
**  o 
g  a. 
•  "<  *» 

.2  §, 
.■ti  *» 
c  S 

l-H  O 


a 

6 

3 

3 


CO 

O 

Ch 

V 

N 

II 

a 

i 

I 

3  tp 
CO  ov 


O 

A 

(-■ 

v 

■•-> 

a 

3 

O 

o 


c4 

6 

5 

§ 

M 

0^ 

0) 

N 

'(3 

6 

M 

O 

2 


•o 

3 

3 

J3 

■u 

2 

.2* 

‘C 

u 

CO 

3 

•s 

o 


'2  S 

C  Ut 
O 

*-•  e 
♦  ^ 

+  3' 

hh 

'O  o* 

(-■  — 
II  II 

a2 

if 

^  GO 


+ 

•«-» 

X 


H ::" 

3'' 2 

0? 

^  a 

%  'B 

pLi 

x-a 

a  CO 

*& 

>«■*•«  ( ) 

■t  * 

CO 

1 

+  2 

X  X 

*9 

' — ’A! 

3 

u 

H  .2* 
c:3  ^ 

W  CO 

1 

II  V 

rand 

ace( 

4> 

4.9 

2 

+9  33 

a  ” 
•c  + 

^  CO 

4>  .s 

II 

••  * 
a^ 

* 

0 

CO  • 

X  • 

‘C  -f- 

^  3 

33  4J 

§  ^ 

J  H 

'b 

49 

SS 

^  3 

2 

a 

3 

4.9  4> 

2  a 

X  1*4 
$-1 

a:;  S 

H  S 

0  a 

•  1^  ^ 

■S  a 

> 


4)  CO 
CO 

J3  * 

o 

'O  'rr' 

(U  2 

-M  ,3 

<2  a 

%l 

CO 

4)  2 

■tj  •*» 

3  a 

S  «-> 
5  “ 
O  II 


JS  ^  ^  ^ 


A-4 


scripts  =  scriptN;  %  initialize  scripts 
for  sindex  =  l:l:k 

%  The  Ptemp  are  complex  vectors  such  that  the  real  part  and  the  complex  part  of 
%  each  of  their  entries  are  normally  distributed  with  means  0.0  and  variances  1.0 


Normalizing  Ptemp  to  obtain  script?  such  that  (scriptP’+scriptP)/T  =  dSNR(sindex) 
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cpower  =  zeros(l,numpts); 
cMtemp  =  scripts  ♦  scriptS’; 


gindex= 1:1:  numpts; 

angle  =  -90  -I-  l*(gindex  -1); 

omega  =  pi  *  sin(angle*pi/180); 
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sumtemp  =  sumtemp  -f-  mpower; 


plot(angle,  10*logl0(cpower),’ 
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plot(angle,  10*logl0(avgpower),’g’) 
hold  off 

%  END  OF  CODE 
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February  19, 1992 


Dr.  C.  V.  Tran 

Naval  Oceans  Systems  Center 
Code  761  (Topside) 

San  Diego,  CA  92152-5000 


Dear  Dr.  Tran: 

With  your  report,  computer  runs,  and  calculations  all  before  me  in  NOSC  Technical 
Document  2226,  October  1991  on  ‘^Beamforming  with  a  Limited  Number  of  Snap¬ 
shots,”  I  can  now  clearly  see  why  the  JASON  beamforming  algorithm  investigated 
in  your  report  appears  to  perform  so  poorly  when  compared  to  conventional  Bartlett 
and  also  MVDR  with  diagonal  loading. 

I  believe  the  distortion  arises  from  a  misinterpretation  of  the  JASON  algorithm,  and 
from  your  choice  thereby  of  the  measure  of  its  effectiveness.  I  will  explain  in  more 
detail  momentarily. 

I  must  also  admit  to  being  disappointed  at  seeing  nothing  in  your  report  about  the 
JASON  work  on  projected  length  of  steering  vectors  into  signal  subspace  -  you  made 
some  computer  runs  on  these  questions  whidi  showed  very  interesting  behavior.  These 
considerations  did  form  an  integral,  though  incomplete,  part  of  the  JASON  algorithm, 
and  might  have  stimulated  further  research.  While  I  personally  intend  none  further, 
I  fear  your  report  in  its  present  form  will  end  all  research  in  these  directions. 

Let  me  begin  with  the  misinterpretation.  You  reinterpreted  the  JASON  algorithm 
as  just  another  beamforming  procedure,  and  appraised  its  performance  in  a  fairly 
standard  fashion.  You  began  by  noting,  as  does  one  of  the  key  arguments  in  the 
JASON  report,  that  if  the  look  direction  steering  vector  7  is  in  the  range  space  of  the 
covariance  matrix  i2,  (and  this  is  the  typical  state  of  affairs  when  signal  is  present  in 
direction  7),  then  the  weight  vector  for  MVDR  beamforming  can  also  be  assun^  to 
lie  in  the  range  space.  Hence  one  can  minimize  w*Rw  subject  to  w*'y  —  1,  with  the 
additional  stipulation  that  tv  is  in  the  range  space  of  R.  This  leads,  as  the  JASON 
report  happily  notes,  to  inverting  a  substantially  smaller  matrix  than  the  conventional 
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diagonally  loaded  MVDR.  But  now  you  go  on  to  use  the  actual  constrained  minimum 
of  w*Rw  with  the  additional  stipulation  above  as  the  array  response  for  any  look 
direction  of  JASON  beamforming.  No  sudi  suggestion  or  clsdm  is  made  in  the  JASON 
report. 

With  this  measure  of  array  response  selected,  the  JASON  algorithm  shows  a  bizarre 
inverted  behavior,  particularly  for  a  small  number  of  snapshots.  While  you  explained 
this  behavior  in  your  report,  I  believe  I  have  a  more  revealing  explanation. 

One  must  ask,  how  does  your  choice  of  array  response  behave  when  the  look  direction 
steering  vector  is  not  in  the  range  space  of  R,  the  state  of  affairs  which  strongly  sug¬ 
gests  no  signal  is  present  in  that  direction.  Well  now,  a  very  simple  calculation  shows 
that  the  response  scales  like  the  reciprocal  of  the  squared  length  of  the  projection  of 
the  steering  vector  into  the  range  space. 

For  the  simulations  you  have  chosen,  with  a  uniform  linear  array,  the  set  of  discrete 
look  direction  steering  vectors  are  almost  perpendicular  to  each  other.  The  direction 
perpendicular  to  the  signal-present  directions  are  thus  only  feebly  present  in  the 
signal  subspace  matrix,  or  equally  in  the  covariance  matrix,  and  so  the  scaling  factor 
is  large  in  the  no-signal  present  directions.  As  the  number  of  snapshots  increases, 
noise  diffuses  into  all  the  look  directions,  and  the  distortion  is  reduced. 

(I  would  not  expect  any  kind  of  beamforming  with  geometrically  dispersed  sensors 
to  work  very  well  if  the  number  of  snapshots  is  much  less  than  twice  the  number  of 
“interesting”  signals  present.  In  your  simulations,  you  are  effectively  getting  super¬ 
resolution  for  Bartlett  end  MVDR  because  the  steering  vectors  are  almost  perpen¬ 
dicular  to  each  other.) 

A  more  reasonable  measure,  then,  of  the  JASON  array  response  would  simply  be 
to  remove  the  scaling  factor.  In  the  notation  of  your  report  this  amounts  simply  to 
multiplying  your  array  response  measure  for  look  direction  7  by  the  factor: 

7*5(5'*5)-^5*7 

7*7 

I  feel  certain  that  the  thus  modified  JASON  array  response  will  be  quite  a  bit  more 
normal  in  appearance.  Yet  I  do  not  like  this  procedure  still,  because  it  is  ad-hoc.  The 
whole  point  of  the  JASON  algorithm  was  to  devise  an  objective  statistical  methodol¬ 
ogy  for  handling  beamforming  with  a  limited  number  of  snapshots.  The  test  we  de¬ 
rived  has  a  measure  of  detect  which  scales  inversely  with  the  length  (not  the  squared 


2 


length)  of  the  projection  of  steering  vector  into  signal  subspace,  and  so  will  not  misbe¬ 
have  quite  as  badly  as  your  array  response  measure.  But  as  your  know,  we  were  still 
very  much  concerned  (with  good  reason)  with  the  large  variance  of  the  detect,  and 
were  trying  to  bring,  objectively,  the  length  of  the  projection  of  the  steering  vector 
into  the  statistical  picture.  We  had  one  sug^tion  for  doing  so,  which  is  far  from 
being  the  end  of  the  matter,  and  it  is  here  that  more  research  needs  to  be  done. 

There  is  another  approach  to  your  measure  of  JASON  array  response,  which  leads 
to  conclusions  similar  to  those  I  have  outlined  above.  If  one  knows  that  7  is  in  the 
ran^c  space  of  i2,  then  as  already  noted,  the  minimum  of  w*Rw  subject  to  10*7  =  1 
may  be  sought  with  w  in  the  range  space  of  R.  But  suppose  that  the  look  direction 
is  not  in  the  range  space.  Then  it  is  easy  to  see  that  the  constrained  minimum  above 
is  simply  zero.  So  one  might  also  have  set  the  JASON  array  response  to  zero  in  these 
directions.  But  we  have  only  a  discrete  imperfect  set  of  look  directions,  and  a  noisy 
estimate  of  R,  so  it  is  natural  to  shade  off  all  responses  according  to  the  length  of 
the  projected  vectors  in  range  space,  perhaps  in  a  manner  similar  to  that  described 
earlier. 

I  hope  you  will  be  able  to  perform  some  additional  computations  for  the  simulations 
described  in  your  report,  which  you  will  circulate  in  a  supplemental  report,  showing 
more  accurately  the  true  potential  of  the  JASON  algorithm,  and  which  may  addi¬ 
tionally  stimulate  further  research.  I  will,  of  course,  be  happy  to  discuss  any  details 
on  further  simulations  with  you. 

I  was  very  pleased  when  NOSC  dedded  some  time  back  to  explore  the  JASON  al¬ 
gorithm,  and  I  am  grateful  for  all  the  work  you  nave  carried  out  to  date.  Perhaps 
I  should  have  paid  closer  attention  along  the  way,  but  as  you  can  see  now,  I  am 
concerned  that  the  current  report  create  a  balanced,  accurate,  and  forward  looking 
description  of  the  JASON  (and  JASON-like)  algorithm.  To  this  end,  I  believe  it  vital 
that  further  detivils  be  circulated. 

Accordingly,  if  you  will  not  be  able  to  supplement  your  origin^  report,  I  am  going  to 
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ask  you  to  distribute  a  copy  of  this  letter  to  all  the  recipients  of  the  original. 
Very  truly  yours, 


cc:  Curtis  G.  Gallan 
Gerald  L.  Mohnkem 
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